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ABSTRACT 


A  method  of  measuring  laser  beam  extinction  in  the 
atmosphere  over  the  ocean  has  been  devised  and  satis¬ 
factorily  utilized.  An  attempt  has  been  made  to  correlate 
the  observed  values  with  existing  meteorological  conditions. 

The  apparatus  consisted  of  equipment,  already  being  utilized 

✓ 

for  scintillation  measurements,  modified  by  the  addition 
of  a  source  intensity  monitor.  Several  calibration  trials 
and  three  over-water  experiments  were  conducted.  A  value 
on  the  order  of  5x10”^  meters  was  found  for  the 
extinction  coefficient  in  clear  weather. 
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I.  PROBLEM  DESCRIPTION 


A.  GENERAL 

The  increasing  utilization  of  laser  technology  in  long- 
range  applications  such  as  communications  links ,  weapons 
designators /range  finders,  and  high  energy  laser  weapons 
makes  it  desirable  to  have  the  capability  of  predicting 
the  effectiveness  of  the  particular  device  at  the  termina¬ 
tion  of  the  path  through  the  atmosphere.  A  highly  collimated 
beam  of  monochromatic  light,  such  as  laser  light,  is  affected 
by  many  phenomena  among  which  scintillation,  beam  wander, 
beam  spread,  and  attenuation  are  only  a  few.  The  develop¬ 
ment  of  a  procedure  by  which  the  latter,  commonly  called 
extinction,  could  be  measured*  in  a  path  through  the  atmos¬ 
phere  directly  above  the  ocean  surface  was  the  goal  of  this 
experiment. 

Atmospheric  attenuation  (extinction)  is  generally 
caused  by  two  phenomena:  absorption  by  the  gaseous  mole¬ 
cules,  and  scattering  by  haze,  fog,  rain,  and  other  aerosols. 
The  absorption  problem  is  minimized  during  system  design 
by  selecting  a  lasing  medium  that  produces  light  of  the 
proper  wavelength  so  as  to  be  within  one  of  the  several 
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transmission  "windows*'  in  the  atmosphere.  Thus,  the  degree 
of  extinction  becomes  a  function  of  the  meteorological 
conditions . 

Ideally,  when  sufficient  data  has  been  obtained,  a 
system  operator  would  only  have  to  observe  the  existing 
meteorological  conditions  and  thus  predict  the  effective¬ 
ness  of  his  equipment.  Extinction  measurements  over 
terrain  have  previously  been  made  /lief.  l7.  However,  such 
measurements  over  the  marine  boundary  are  virtually  non¬ 
existent  . 

B.  MATHEMATICAL  DESCRIPTION 

If  a  beam  of  monochromatic  flux  enters  an  attenuating 
layer,  the  Bouger-Lambert  Law  /Ref.  X]  states  that  the 
fliix  at  any  point  within  the  layer  is  given  by 

P=  PoS'^’^  (1) 

where  P=  Flux  remaining  after  the  beam 
has  traversed  a  path  length  %  , 

Po  =  Flux  entering  the  attenuating 
layer, 

=  Extinction  Coefficient, 

•pc  =  Path  length  in  the  attenuating 
layer. 

The  flux  described  in  equation  (1)  above  has  the  units 
of  energy  per  unit  time  over  a  given  area  and  will  be 
further  described. 
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Considering  both  absorption  and  scattering  the  extinc¬ 
tion  coefficient 

becomes  =  o<  ^  (2) 

where  <=*<  =  absorption  coefficient, 

and  ^  =  scattering  coefficient. 

In  this  experiment,  0.6328  micrometer  wavelength  laser 
light  was  utilized  and  since  this  is  within  one  of  the 
transmission  "windows"  no  attempt  was  made  to  distinguish 
the  relative  values  of  *=<  and  ^  ,  but  only  to  determine 
their  combined  value. 

The  relationship  (Equation  1)  holds  for  a  spherical 
wave.  The  detailed  mathematical  development,  following 
references  3  and  4,  presented  in  Appendix  A,  shows  that  a 
laser  source  producing  pure  monochromatic  radiation  can  be 
considered  to  be  a  point  source  of  spherical  waves,  pro¬ 
vided  the  observations  are  made  in  the  paraxial  region,  at 
a  distance  from  the  source  much  greater  than  the  off-axis 
position. 

Thus,  considering  again  equation  (1),  and  assuming  that 
the  effective  source  of  the  radiation  is  at  the  Gaussian 
waist  of  the  diverged  beam  and  radiates  P®  watts  into  a 
solid  angle  ,  then  the  intensity  at  a  position  near  the 

beam  axis  at  a  distance  R  is  given  by 


11 


I 


Po  e 


- 


(3) 


2 

where  I  =  Intensity  at  R  in  Watts /meter 

2 

_ri_  R  =  Area  of  beam  at  distance  R. 
Consider  propagation  in  a  controlled  atmosphere  where 
scattering  aerosols  are  at  a  minimum,  so  that  , 

and  a  calibration  factor  can  be  determined  from 

T  =  (4) 

r2 

where  the  subscript  indicates  all  measurements  taken  during 
the  controlled  trial.  Then  from  (4) 


(5) 


so 
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^  loRt  -e  e 
r2 


(6) 


Taking  the  Natural  Logarithm  results  in 

In  (IR^)  =  In  (I^Ro  ^  -  yU,R.  (7) 

2 

A  plot  of  In  (IR  )  versus  R  will  yield  the  ordinate 
intercept 

In  (I^R^  =  In  C.  (8) 
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So  in  general  for  longer  ranges  where  the  determination 
of  jul,  is  desired 


I 


 Ce 


(9) 


and 


from  which: 


In  IR^  =  In  C  y<.R 


In  C  -  In  IR^. 


(10) 

(11) 


Values  of  In  /iR^J  at  each  value  of  R  can  thus  yield 
a  calculated  value  of  . 

Alternatively,  it  is  only  necessary  to  plot  In  £  IR^^ 
versus  R  for  a  series  of  measurements  at  increasing  ranges 
and  from  this  determine  yU.  as  the  absolute  value  of  the 
slope.  This  assumes,  of  course,  that yu.  remains  constant 
over  the  time  and  distance  of  the  measurements. 

To  obtain  the  necessary  data  for  the  plot  of  In  /  IR^J 
versus  R  the  procedure  briefly  described  below  is  followed: 

A  horizontally  fanned  beam  of  constant  intensity  sweeps 
vertically  across  a  detector  located  at  a  known  distance  R, 
4000  times  each  second.  The  peak  value  of  each  detected 
pulse  is  recorded  and  the  pulse  height  distribution 
(spread  out  because  of  scintillation)  is  obtained.  The 
average  value  of  the  pulse  height  is  calculated  and  taken 
to  be  a  measure  of  the  average  intensity  I,  at  the  distance 
R. 
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II.  EXPERIMENTAL  APPARATUS 


A.  GENERAL 

Since  the  ultimate  goal  of  the  experiment  was  to  measure 
the  extinction  coefficient  over  the  marine  boundary  in  con¬ 
junction  with  simultaneous  meteorological  measurements,  the 
use  of  the  Research  Vessel  ACANIA  was  necessary.  Because 
of  the  limited  availability  of  the  ship's  services,  it  was 
deemed  desirable  to  test  equipment  and  procedures  thoroughly 
prior  to  actual  shipboard  measurements.  This  was  accomplished 
first  in  the  120  meter  long  basement  hallway  of  the  Naval 
Postgraduate  School's,  Spanagel  Hall.  To  accomplish  this 
the  laser  source  and  scanning  equipment  (to  be  described 
in  detail  in  the  next  section)  were  placed  at  one  end  of 
the  hallway,  and  the  detector  and  processing  equipment  were 
moved  the  length  of  the  hallway  in  twenty  meter  increments. 

When  operation  of  the  equipment  in  this  controlled 
environment  was  deemed  satisfactory  the  apparatus  was  re¬ 
configured.  The  laser  source  and  scanning  equipment  were 
mounted  on  a  training  and  elevating  mechanism  in  a  two¬ 
wheeled  trailer.  The  detecting  and  processing  equipment 
was  mounted  in  an  instrumented  bus  (Fig.  1)  which  provided 
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Naval  Postgraduate  School  Mobile  Optical 
Research  Laboratory. 


Interior, Showing  18'lnch  Telescope. 


FIGURE  1. 


INSTRUMENTED  BUS  and 


IB'Inch  Telescope 
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portability,  shelter,  and  110  V.A.C.  power  from  a  3.5  K.W. 
gasoline  engine  driven  generator.  This  configuration  pro¬ 
vided  the  necessary  support  for  extending  the  path  length 
to  300  meters  in  100  meter  increments  along  O'HARE  Avenue 
at  the  Naval  Postgraduate  School  Annex.  Operation  in  this 
manner  was  the  final  preliminary  check  prior  to  shipboard 
operation.  For  shipboard  measurements  the  apparatus  was 
again  re-configured.  The  laser  source  and  scanning  equip¬ 
ment  were  mounted  in  the  bus  and  the  detecting  and  processing 
equipment  were  aboard  ship.  The  reasoning  behind  this  final 
re-configuration  will  become  apparent  in  the  following 
paragraphs . 

B.  DETAILS  OF  EXTINCTION  EQUIPMENT 

The  laser  source,  scanning  equipment,  detection,  and 
processing  equipment  utilized  for  extinction  measurements 
were  also  utilized  for  scintillation  and  aperture  averaging 
measurements;  however,  only  that  equipment  in  direct  support 
of  the  extinction  experiment  will  be  described  in  detail. 

1.  Bus  Mounted  Equipment 

Equipment  described  in  the  following  paragraphs 
was  mounted  aboard  the  instrumented  bus  utilized  by  the 
optical  propagation  group.  The  bus  itself  was  a  converted 
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General  Motors  Corporation  passenger  bus.  The  forward 
one-third  of  the  vehicle  contained  storage  cabinets  and 
workbench  space  and  the  drivers  compartment.  The  after 
two -thirds  of  the  vehicle  contained  electronic  equipment 
racks,  workbench,  and  an  18-inch  telescope  mounted  on  a 
platform  capable  of  being  raised  clear  of  the  floor  on 
three  hydraulic  camper  jacks  extending  through  the  floor 
to  the  ground.  This  arrangement  provided  a  stable  mount 
for  the  telescope  and  all  equipment  mounted  on  it  and  was 
thus  isolated  from  any  vibration  caused  by  personnel  moving 
about  in  the  vehicle.  A  portable  3.5  kilowatt  generator 
provided  110  V.A.C.  power  to  all  equipment, 
a.  Laser  Source 

The  laser  utilized  as  a  source  for  the  extinc¬ 
tion  measurements  was  a  0.6328  micrometer  He-Ne  laser 
manufactured  by  C.W.  Radiation  INC.  The  output  of  the  laser 
was  fitted  with  a  diverging  lens  controlled  by  a  micrometer 
screw  adjustment.  The  laser,  scanning  equipment,  and  source 
intensity  monitor  were  mounted  on  a  common  plate  attached 
to  the  right -hand- side  of  the  bus  mounted  18-inch  telescope 
(Fig.  1) .  This  provided  a  means  by  which  the  laser  could 
be  trained,  elevated,  and  aimed  at  the  shipboard  mounted 
detectors . 
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b.  Scanning  Equipment 

(1)  Physical  Description.  Mounted  directly 
in  front  of  the  laser  was  a  mirror  positioned  at  a  45° 
downward  angle  (Fig.  2) .  This  mirror  performed  two  func¬ 
tions  .  First,  it  directed  the  beam  down  onto  the  scanning 
mirror.  Secondly,  it  was  capable  of  being  bent,  by  means 
of  a  micrometer  screw-adjustable  knife-edge  placed  in  the 
vertical  axis,  so  as  to  produce  a  "FANNED”  beam  elongated 
in  the  horizontal  plane.  This  arrangement  was  capable  of 
producing  a  beam  at  the  target  ship  sufficiently  wide  so  as 
to  provide  uniform  coverage  of  the  detector  as  the  ship 
moved  and  swung  about  its  anchor. 

The  actual  scanner  consisted  of  a  small, 
1.0  cm  X  1.0  cm,  mirror  mounted  at  a  45°  upward  angle 
^ef.  6,  p.  4Q7  .  This  mirror,  attached  to  two  piano  wires 
and  mounted  between  the  poles  of  a  strong  magnet,  was  set 
into  oscillation  and  driven  by  a  power  amplifier.  The 
amplifier  controls  were  adjustable  so  as  to  provide  a 
maximum  vertical  scan  angle  of  14  milliradians  at  a  fre¬ 
quency  of  approximately  2000  HZ.  The  result  of  this  beam 
processing  was  to  allow  detection  of  maximum  intensity 
independent  of  the  ship's  vertical  motion. 
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Cylindrical  Mirror 
(micrometer  adjustment  provides 


19 


Retractable  Mirror 

^  _  ,  (shown  in  position  for 

FIGURE  2  Diagram  of  Beam  Path  sampling  beam  intensity) 


(2)  Mathematical  Description.  The  following 


is  a  brief  analysis  of  the  scanned  beam  in  terms  of  the 
vertical  angular  width  of  the  beam  and  the  scanning  rate 
of  the  oscillating  mirror.  This  permitted  a  prediction  of 
the  time  width  of  the  detected  pulse  which  was  then  found 
to  be  in  agreement  with  the  measured  value.  To  determine 
this  the  beam  was  fanned  and  both  the  unscanned  and  scanned 
beam  sizes  were  measured  at  a  convenient  distance  from  the 
scanning  mirror.  Refer  to  Figure  3  for  the  details  of  the 
geometry;  where  ^  (the  distance  from  the  scan  mirror)  = 
2316.5  cm,  h^  (the  unscanned  beam  height  at  /  )  =5.5  cm, 
and  h.2  (the  scanned  beam  height  at  /  )  =  26  cm. 


From  this  c<  ,  the  unscanned  beam  angle 


of  coverage,  equals  2.37  milliradians  and  2-0^,  the  angle  of 
coverage  of  the  scanned  beam,  equals  11.22  milliradians. 
Thus,  -0^  =5.61  milliradians. 


Now  the  pulse  width  is  given  by 


where 


/V  -  ^  ^ 


(12) 


scanning  mirror  by  the 


detector. 


and  ~  maximxim  angular  scan  rate. 
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FIGURE  3  Geometry  of  Pulse  Width  Determination 
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Since 


for  0  where 


/"  =  f 


D  equals  the  height  of  the  detector. 


Thus , 


r- 


c<  + 


(13) 


For  large  values  of  ^  ,  the  second  term 
in  the  numerator  of  (13)  becomes  negligible  thus  rendering 
the  time  width  of  the  pulse  nearly  constant  with  range. 


Now 


-e-  ^ 


5;o  CD  t) 


or 


d-e 


cj>  -Q-  Cos  cD)t: . 


(14) 


(15) 


For  maximum  value  of  -O-  ,  Cos  uj't  ~  dL 
and 


(16) 


with 


^  3.000  ) 

and  ■%  =5.61  milliradians ,  as  before; 

Thus,  =  70.5  radians /sec. 

For  a  detector  of  two  centimeter  diameter 
and  Jl  as  before,  equation  (13)  yields  a  time  width  of  the 
pulse,  ^  ,  of  45.8  microseconds .  This  agreed  with 

oscilloscope  measurements  of  the  detected  pulses. 
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Figure  4  shows  the  scanning  equipment  and 
mounting  details;  the  entire  assembly  was  referred  to  as 
the  "FAN-SCAN," 

c.  Source  Intensity  Monitor 

Periodically  throughout  the  extinction  measure¬ 
ments  it  was  necessary  to  sample  the  laser  output  as  the 
beam  left  the  scanning  mirror  to  ensure  that  no  intensity 
degradation  was  present  at  the  source.  To  facilitate  this, 
the  Source  Intensity  Monitor  was  devised  (Fig.  4).  An 
adjustable  mirror  was  installed  which  could  be  swung  into 
position  to  interrupt  the  outgoing  beam  and  direct  it  back 
to  the  monitor  detector  mounted  above  and  behind  the  source 
laser.  The  detector,  a  PIN-5D  Photodiode  was  mounted  in 
an  adjustable  mount  to  permit  fine  alignment  of  the  sensi¬ 
tive  surface  with  the  sampled  beam.  The  output  of  the 
detector  was  displayed  on  a  Monsanto  OS  226  P,  oscilloscope. 

2 .  Shipboard  Mounted  Equipment 

The  detection  and  signal  processing  equipment  aboard 
the  R.V.  ACANIA  included,  detector,  gyro-stabilized  detector 
mount,  PAR  113  Low-Noise  Differential  Preamplifier,  Demodu¬ 
lator,  Hewlett-Packard  7562A  Log  Voltmeter,  Victoreen  PIP- 
400  Multichannel  Analyzer,  and  Hewlett-Packard  HP  9810 
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Laser  Source  and  Scanning  Equipment  Mounted 
on  ISHnch  Telescope. 


CloseMip  Showing,  from  left, Monitor  Detector, 
Laser, Cyl indrical  Mirror,  Oscillating  Mirror, 
and  Retractable  Mirror. 


FIGURE  A.'FAN-'SCAN  and  Source  Intensity  Monitor  Arrangement 
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computer  and  associated  plotter.  The  following  paragraphs 
describe  in  detail  the  arrangement  and  fxinction  of  each  of 
the  above. 

a.  Detector 

The  detector  utilized  for  the  extinction  measure¬ 
ments  was  a  TIX  L69  Silicon  Avalanche  Photodiode  reversed 
biased  to  159  volts.  This  provided  an  acceptable  Signal- 
To-Noise  Ratio  at  the  extended  ranges.  The  detector  was 
mounted  in  a  holder  fitted  with  lens  and  an  aperture  adjust¬ 
able  to  4.2  cm  diameter.  For  more  detailed  discussion  of 
the  detector  see  ZRef.  77. 

b.  Gyro  Stabilizer 

The  detector  was  mounted  on  the  gyro  stabilized 
mount.  The  stabilization  was  necessary  to  maintain  align¬ 
ment  of  the  detector  with  the  laser  source  ashore. 

The  gyro  was  an  inertially  stabilized  platform 
utilizing  two  iron  flywheels  mounted  on  each  end  of  a  1/2  hp 
capacitor  start  electric  motor  (Fig,  5).  The  shore  based 
laser  source  utilized  for  the  extinction  measurements  also 
served  as  the  beacon  onto  which  the  gyro  feedback  circuit 
locked.  Feedback  was  provided  by  a  four-quadrant  position 
and  intensity  photodetector  with  associated  electronics. 
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Gyro  unit;Showing  top  mounted  Avalanche 
photodiode  detector. 


Signal  processing  equipment;  PARIIS’sCleftl; 
equipment  in  the  rack  ( top  to  bottom)  ;  Demodula- 
tor,  Log  Voltmeter/ Amplifier,  Oscilloscope, and 
Multi-channel  Analyzer. 


Hgure  5.  Signal 


Detection  and  Processing  Equipment 
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A  small  aiming  telescope  and  "JOY  STICK"  control 
were  mounted  on  the  left  side  of  the  gyro  to  aid  the  gyro 
operator  in  aiming  the  control  detector  at  the  laser  source. 
The  gyro  platform  was  positioned  on  the  upper  deck  of  ACANIA 
and  was  sheltered  from  the  elements  by  an  aluminum  tool 
shed  open  on  one  side.  The  signal  processing  equipment 
described  below  was  bench  mounted  in  the  "DRY  LAB"  on  the 
main  deck.  (Fig.  5) . 

c.  Preamplifier 

The  output  of  the  detector  was  fed  through  an 
F.E.T.  follower  (which  served  the  function  of  impedance 
matching)  to  the  PAR  113  Low-Noise  Preamplifier.  This 
provided  high  gain,  10  to  10^,  low  noise  amplification,  of 
the  received  signal.  Adjustable  High  and  Low  Frequency 
Rolloffs  allowed  selection  of  the  bandwidth. 

d.  Demodulator 

The  modulated  signal  was  next  passed  to  the 
Demodulator.  This  device  sensed  the  scanning  beam  pulse, 
followed  the  pulse  to  its  maximum  and  held  that  instanta¬ 
neous  signal  level  in  a  flat-topped  pulse  until  the  next 
maximimi  scanned  across  the  detector  at  which  time  it  held 
that  maximum  level,  and  so  on.  The  output  of  this  circuit 
was  thus  a  series  of  step  functions. 
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e.  Log  Voltmeter 

The  Hewlett-Packard  7562-A  Log  Voltmeter  was 
the  next  instrxmient  in  the  processing  step.  At  this  point 
the  Logarithm  of  the  signal  was  obtained. 

f.  Multichannel  Analyzer 

The  Victoreen,  PIP-400,  Multichannel  Analyzer 
collected,  stored  for  further  readout,  and  displayed  the 
sampled  Log  Voltage  Pulses.  Each  of  the  available  200 
channels  corresponded  to  a  different  pulse  height  and  the 
ntjmber  of  pulses  of  any  particular  height  received  during 
the  sampled  period  were  stored  in  the  respective  channels. 

g.  Computer  and  Plotter 

The  output  of  the  Multichannel  Analyzer  was 
fed  into  a  Hewlett-Packard  9810A  calculator.  The  calculator 
provided  controlling  functions  to  the  associated  plotter. 

The  calculator  sampled  each  channel  of  the  Multichannel 
Analyzer,  stored  the  information  and  plotted  the  pulse 
height  distribution.  The  calculator  then  computed  a  least- 
squares  polynomial  fit  of  a  Gaussian  distributed  curve  to 
the  data  points.  Using  the  position  of  this  curve  with 
respect  to  read-in  calibration  points  and  the  Standard 
Deviation  of  the  curve,  the  average  intensity  of  the  re¬ 
ceived  signal  was  computed.  See  Appendix  B  for  a  detailed 
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analysis  of  the  procedure  for  obtaining  the  average  pulse 
height  from  the  log  voltage  distribution.  For  a  detailed 
description  of  the  computer  program  see  Ref.  8.  Figure  6 
is  a  sample  printout  from  the  system. 

h.  HTjmidity  Measuring  Equipment 

The  Hygrodynamics  Digital  I  hygrometer  indicator 
and  the  Dunmore-t3rpe  Lithium  Chloride  sensor  were  used  to 
measure  the  Relative  Humidity.  The  equipment  operates  on 
the  basis  of  resistance  change  within  an  electrolytic 
solution  generating  a  reference  voltage  variance  which  is 
proportional  to  the  relative  hvimidity  change.  For  further 
details  see  Ref.  9. 
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FIGURE  6  Sample  Computer  Output  Showing  "Best  Fit"  Gaussian  Curve 
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III.  THE  EXPERIMENT 

A.  EXPERIMENTAL  GOALS 

The  original  goals  of  the  experiment  were  threefold. 
First,  the  main  goal  was  the  development  of  the  equipment 
and  procedures  to  measure  extinction  and  to  demonstrate 
that  extinction  could  be  measured  in  the  ocean  environment. 
Secondly,  it  was  desirable  to  correlate  the  observed  ex¬ 
tinction  with  the  measured  quantities  of  airborne  particu¬ 
late  matter  and  water  vapor  measured  as  humidity.  Thirdly, 
it  was  desired  to  measure  the  extinction  caused  by  fog  and 
rain.  This  last  goal  was  not  achieved  since,  during  the 
scheduled  shipboard  trials,  the  weather  was  persistently 
clear. 

B.  EXPERIMENTAL  PROCEDURE 

1.  22-23  January  1976 

On  these  dates  the  first  equipment  trials  were 
performed  in  the  basement  hallway  of  Spanagel  Hall,  Naval 
Postgraduate  School.  The  laser  source  and  monitoring 
equipment  were  mounted  on  a  portable  equipment  rack  and 
positioned  at  distances  of  15.8  m,  31.0  m,  and  46.3  m  from 
the  receiving  equipment  which  was  kept  stationary.  The 
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purpose  of  these  trials  was  to  develop  technique  and 
initial  equipment  settings.  At  that  time  the  signal 
detector  was  a  PIN  Photodiode  which,  for  reasons  to  be 
explained  in  Section  III.B.5,  was  replaced  by  the  Avalanche 
Photodiode  previously  described.  In  addition  to  the  devel¬ 
opment  of  techniques  and  equipment  settings,  these  trials 
established  the  Laser  Source  Strength  Monitor  reference 
reading  at  31  mv. 

2.  26  January  1976 

On  this  date  the  previous  trials  were  repeated  with 
similar  results.  In  addition,  the  effects  of  "Fan-Scan" 
amplifier  settings  and  PAR  amplifier  settings  were  measured. 
It  was  determined  that  "Fan-Scan"  amplitude  had  a  definite 
effect  on  the  monitor  pulse  height  and  an  amplifier  setting 
of  1.2  was  chosen  as  the  best  setting  since  the  response 
at  settings  surrounding  this  point  was  relatively  constant. 
The  PAR  settings  arrived  at  were:  (1)  Low  Frequency  Rolloff 
100  HZ,  (2)  High  Frequency  Rolloff  300  KHZ.  These  settings 
provided  the  cleanest  signal  response  with  minimijm  pulse 
shape  distortion. 

3.  2  February  1976 

On  this  date,  the  first  concrete  attempt  was  made 
to  determine  the  equipment  response  to  transmission  through 
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an  atmosphere  containing  essentially  no  water  droplets. 

The  trial  was  again  made  in  the  basement  hallway  at  trans¬ 
mission  distances  of  15.8  m,  31.0  m,  46.2  m,  61.5  m,  92.0  m 
and  122.5  m.  All  equipment  settings  were  maintained  as 
previously  described;  however,  in  this  trial  the  laser 
source  was  left  stationary  and  the  detector  and  signal 
processing  equipment  were  positioned  at  the  prescribed 
ranges.  This  procedure  was  utilized  since  it  facilitated 
more  rapid  detector  alignment.  The  results  of  this  trial 


are  plotted  in  Figure  7.  The  observed  value  of  In 


determined  here  was  to  be  utilized  as  In 


in  the  cal¬ 


culation  of  the  extinction  coefficient  as  a  result  of  the 
data  to  be  taken  during  the  shipboard  trials  of  the  week 
of  29  March  1976. 

4.  18  March  1976 


It  was  decided  that,  prior  to  the  29  March  Acania 


trip,  it  would  be  advantageous  to  have  some  measurements 
at  longer  distances  than  could  be  achieved  in  the  basement 
hallway.  For  the  calibration  trial  on  this  date  the  laser 
source  and  monitoring  equipment  was  mounted  in  the  small 
two-wheeled  trailer,  and  positioned  at  the  western  end  of 
O' Hare  Avenue  (NAVY  ANNEX)  in  close  proximity  to  the  fuel 
station  so  that  115  volt  A.C.  power  could  be  readily 
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Spanagel  Hall  Basement 
2  February  1976 


% 


\ - H - 1 - V 


m 


c^yn  ui 


H - V 

CM  ^ 


34 


100  200 
RANGE  IN  METERS 

FIGURE  7  Graph  of  In(IR^)  vs.  R  for  2  February  1976 


obtained  from  the  distribution  panel  located  there.  (See 
Figure  8) .  The  detector  and  signal  processing  equipment 
was  located  in  the  instrumentation  bus,  which  was  moved 
and  positioned  at  ranges  of  91.0  m,  206.0  m,  and  336.0  m 
from  the  source.  (See  also  Figure  8).  The  results  of  this 
trial  are  plotted  in  Figure  9.  These  results  together  with 
the  results  of  the  2  February  trials  were  to  be  utilized 
with  the  data  collected  aboard  Acania  during  the  week  of 
29  March. 

5 .  Week  Beginning  29  March  1976 

During  this  week,  several  extinction  measurements 
were  recorded  aboard  R.  V.  Acania.  The  majority  were  taken 
with  the  ship  at  anchor  at  a  range  of  approximately  0.8 
nautical  miles  off  Point  Pinos.  Each  day  the  ship  would 
arrive  on  station  and  anchor  using  radar  ranges.  When  the 
ship  settled  and  rode  to  the  anchor,  the  exact  range  from 
bus  to  ship  was  determined  by  triangulation  using  surveyor's 
transits.  It  was  during  this  trial,  on  30  March  1976,  when 
an  attempt  was  made  to  obtain  extinction  measurements  at 
consistently  increasing  ranges,  that  it  was  determined  that 
to  do  so  would  require  the  use  of  the  Avalanche  Detector; 
since,  at  ranges  greater  than  two  kilometers,  the  Signal- 
to-Noise  Ratio  obtained  with  the  PIN  Photodiode  was 
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FIGURE  8  Drawing  of  Experimental  Arrangement  Naval  Postgraduate  School  Annex 


I 


FIGURE  9  Graph  of  In(IR^)  vs.  R  for  18  March  1976 


insufficient.  From  this  determination,  it  was  decided  that 
a  "calibration"  trial,  similar  to  the  18  March  trial,  would 
have  to  be  made  using  the  Avalanche  Photodiode.  Since  this 
detector  produces  an  extremely  large  signal  resulting  in 
saturation  at  close  ranges,  it  was  decided  to  again  work 
at  the  Navy  Annex  site  rather  than  the  Spanagel  Hall  basement. 

6.  19-20  April  1976 

On  these  dates  the  "calibration"  trials  were  con¬ 
ducted  utilizing  the  Avalanche  Photodiode.  The  equipment 
was  again  configured  as  described  previously  in  Paragraph  4. 
The  source  strength  monitor  at  this  time  also  produced  a 
reading  of  31  iitv  indicating  that  the  laser  output  was  re¬ 
maining  constant.  The  results  of  the  trials  described  here 
are  plotted  in  Figure  10.  Again,  the  observed  value  of 
In  I  determined  here  was  to  be  utilized  as  In  |cj 

in  conjunction  with  the  data  ,to  be  collected  during  the 
final  shipboard  trials  to  be  conducted  the  following  week. 
Also,  during  this  calibration  trial,  anticipating  the  need 
to  use  a  large  aperture  opening  at  the  longer  ranges,  a 
trial  was  conducted  to  determine  the  effect  on  the  response 
of  the  detector  as  a  function  of  aperture  size.  This  was 
conducted  at  the  longest  available  range  (336.0  m)  since 
opening  the  aperture  greater  than  0.4  cm  diameter  at  the 
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FIGURE  10  Graph  of  InClR^)  vs.  R  for  19-20  April  1976 


closer  ranges  resulted  in  saturation.  The  results  of  this 
trial  are  plotted  in  Figure  11.  As  indicated  the  response 
is  linear  to  an  aperture  radius  of  1.1  cm,  at  which  point 
saturation  set  in. 

7.  28-29  April  1976 

On  these  dates,  the  final  data  collection  reported 
on  in  this  thesis  was  completed.  Simultaneous  meteorological 
data  was  collected  at  the  detector  site  aboard  ship.  The 
laser  source  was  mounted  on  the  18  inch  telescope  in  the 
bus.  The  bus  was  positioned  at  the  water's  edge  on  Point 
Pinos.  The  ship  was  positioned  north  of  this  location  at 
various  ranges  from  839  meters  to  6748  meters.  Three 
factors  arose  that  affected  the  effort  to  duplicate  the 
exact  equipment  conditions  of  the  calibration  trials.  First, 
because  of  the  heavy  electrical  load  on  the  3.5  kw  generator 
resulting  in  a  low  line  voltage,  the  "Fan-Scan"  amplifier 
setting  had  to  be  increased  to  2.5  to  maintain  oscillation 
of  the  scanning  mirror.  Secondly,  the  source  intensity 
monitor  indicated  some  laser  degradation  as  a  result  of 
the  low  voltage  condition.  The  maximum  monitor  reading 
was  22.0  mv  when  the  trial  first  began,  then  within  30 
minutes  of  operation  the  monitor  dropped  to  18.5  mv  and 
remained  there  throughout  the  remainder  of  the  trials. 
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FIGURE  11  Graph  of  Average  Intensity  vs.  Aperture  Size 


This  was  down  from  a  constant  31.0  mv  obtained  during  the 
calibration  trials.  Thirdly,  because  of  the  long  ranges 
involved,  it  was  decided  that  the  best  procedure  would  be 
to  set  the  aperture  size  at  4.0  cm  diameter  and  keep  that 
constant  setting  throughout  all  remaining  trials. 

The  equipment  settings  were  as  follows: 

PAR  Amplifier:  Low  Frequency  Roll  Off,  100  Hz., 

High  Frequency  Roll  Off,  300  kHz; 

"Fan-Scan":  Amplifier,  2.5, 

Cylindrical  mirror  micrometer,  7; 
Laser  divergence  micrometer,  5; 
Detector  Aperture,  4.0  cm  diameter. 

The  general  experimental  procedure  for  these  two 
days  was  to  maneuver  the  ship  into  position  north  of  the 
bus  at  the  desired  range.  When  on  station  the  gyro  stabi¬ 
lized  detector  would  be  aimed  at  the  laser  source.  When  a 
usable  signal  was  obtained  the  data  collection  would 
commence.  The  Multichannel  Analyzer  would  be  allowed  to 
collect  a  sufficient  number  of  counts  so  as  to  produce  a 

I 

nearly  Gaussian  distribution.  This  usually  took  from  one 
to  three  minutes  depending  upon  distance  from  the  source. 
Since  it  was  nearly  impossible  to  hold  the  ship  at  the 
same  range  throughout  data  collection,  approximately  half 
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way  through  the  collection  process  the  ship's  range  from 
the  source  would  be  determined  by  triangulation  with  the 
surveyor's  transits.  As  soon  as  the  data  collection  was 
completed  the  ship  would  transit  to  the  next  desired 
position  and  the  process  would  be  repeated.  While  the  ship 
was  moving  to  the  new  position  the  data  was  fed  into  the 
HP  9810  Computer,  processed,  and  plotted  on  the  associated 
plotter;  thus  freeing  the  Multichannel  Analyzer  to  collect 
data  at  the  next  position.  At  the  closer  ranges  the  gyro 
would  "lock-on"  the  laser  signal;  however,  at  ranges  greater 
than  about  3000  meters  the  signal  strength  was  not  sufficient 
for  "lock-on"  thus  the  gyro  operator  controlled  the  aiming 
of  the  detector  manually  with  the  "Joy-Stick"  control.  The 
results  of  the  data  collected  and  the  correction  factors 
resulting  from  the  equipment  changes  previously  described 
are  reported  in  Section  IV,  RESULTS  AND  CONCLUSIONS. 
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IV.  RESULTS  AND  CONCLUSIONS 


A.  GENERAL 

Since  the  major  goal  of  the  experiment,  as  stated 
earlier,  was  the  determination  of  the  extinction  coefficient 
within  the  marine  boundary  layer,  the  actual  data  gathering 
revolved  about  the  schedule  of  the  Research  Vessel  Acania. 
Because  of  this  and  the  discovery  on  30  March  of  the  long 
range  limitations  of  the  PIN  Diode  detector,  the  data 
breaks  down  logically  into  two  distinct  experimental  evolu¬ 
tions:  (1)  the  30  March  experiment  and  trials  leading  up 

to  it,  and  (2)  the  28-29  April  experiment  and  the  calibra¬ 
tion  trials  leading  to  it. 

To  facilitate  the  determination  of  the  extinction  co¬ 
efficient  ,  >  considering  the  paucity  of  long-range  data, 

the  Natural  Logarithm  of  the  quantity  (IR^)  was  plotted 
versus  the  range,  R,  for  each  data  point  obtained.  A 
straight  line  was  drawn  through  these  data  points  utilizing 
the  ordinate  intercept  and  slope  as  determined  from  the 
Principle  of  Least  Squares  curve  fitting  /Ref.  1^.  The 
absolute  value  of  the  slope  of  this  curve  is  equal  to » 
the  extinction  coefficient. 
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The  calibration  trials  served  to  provide  short-range 
data  and  thus  determine,  to  a  large  degree,  the  ordinate 
intercept  of  the  plot;  whereas,  the  long-range  data  obtained 
in  the  desired  ocean  environment,  from  Point  Pinos  to  Acania, 
provided  the  information  for  the  slope  and  consequently • 

Prior  to  plotting  the  data,  the  following  corrections 
were  made  where  necessary: 

1.  Source  Intensity  Correction 

Since  the  majority  of  the  data  was  taken  with  a 
laser  source  intensity  monitor  reading  of  31.0  mv,  the 
pulse  height  data  taken  while  the  laser  intensity  monitor 
reading  differed  from  this  was  corrected  either  up  or  down 
according  to  the  ratio  of  the  two  readings. 

2 .  Aperture  Size  Correction 

As  mentioned  earlier,  it  was  necessary  to  utilize 
a  small  aperture  at  short  ranges  and  a  larger  aperture  at 
long  ranges.  Therefore,  the  data  collected  at  the  longer 
ranges  was  reduced  by  the  ratio  of  the  squares  of  the 
aperture  radii. 

B.  THE  30  MARCH  EXPERIMENT 

The  30  March  1976  experiment  and  the  calibration  trials 
proceeding  it  served  to  establish  the  feasibility  of  making 
measurements  at  increasing  ranges  as  well  as  providing  the 
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initial  measurement  of  the  extinction  coefficient  in  the 
marine  environment.  The  data  is  tabulated  in  Appendix  C. 
Figure  12  is  the  plotted  data  including  the  least  squares 
fitted  straight  line.  The  weather  was  clear  and  the  simul- 

j 

taneous  meteorological  observations  provided  an  average 

i 

relative  humidity  reading  of  0.911  over  the  time  of  the 
shipboard  measurements .  was  6.2  x  10"^  meters 

C.  THE  28-29  APRIL  EXPERIMENT 

The  experiments  conducted  on  28  and  29  April  1976, 
together  with  the  preceeding  calibration  trials,  succeeded 
in  providing  data  at  even  longer  ranges .  The  data  taken 
on  these  two  dates  was  considered  together  for  calculation 
purposes  because  of  the  similar  weather  conditions  and 
equipment  operation.  The  data  is  tabulated  in  Appendix  C. 
Figure  13  is  the  plotted  data  including  the  least  squares 
fitted  straight  line.  The  weather  was  again  clear  and  the 
relative  hxmiidity  on  28  April  was  0.776,  while  on  29  April 
it  was  0.774  giving  an  average  relative  humidity  reading 

of  0.775  over  the  time  of  the  measurements . was  3.8  x 

-4-  -1 

10  meters 
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Pt.  Pinos  to  ACANIA 
30  March  1976 
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FIGURE  12  Graph  of  InClR"^)  vs.  R  for  30  March  1976 
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FIGURE  13  Graph  of  ln(IR2)  vs.  R  for  28-29  April  1976 


D.  CONCLUSIONS 


It  seems  possible  to  draw  three  main  conclusions  from 
these  results. 

First,  it  is  experimentally  feasible  to  measure  laser 
beam  extinction  at  sea  over  the  ocean-air  interface. 

Secondly,  the  extinction  coefficient  for  0.6328  micro¬ 
meter  laser  light  in  clear  weather,  through  the  marine 
environment  is  of  the  order  of  5  x  10”^  meters 

Thirdly,  humidity  may  be  considered  as  a  possible  para¬ 
meter  with  which  to  predict  laser  beam  extinction  in  clear 
weather . 

The  latter  two  items,  of  course,  require  further 
investigation. 
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V.  DISCUSSION  OF  ERRORS 


The  considerations  and  assumptions  discussed  below  were 
utilized  in  an  attempt  to  understand  how  the  uncertainties 
in  the  measured  quantities  affect  the  uncertainty  in  the 
determination  of^^^  . 


(17) 


the  more  useful  form 


Jn  =  Z  c  -yun, 


(18) 


was  used.  Then yuc  was  determined  to  be  the  absolute  value 
of  the  slope  of  the  plot  of  In  IR'^  versus  R.  From  (18) 


Me  -  J^r  -  i/.R 


(19) 


From  the  theory  of  propagation  of  errors  /Ref.  10_7 ; 


From  equation  (19), 


ifi 

c)  C- 


-  _L  ^  J- 

'  -Re.  >  jx  7?X 


and 


_  _  7  A  I  .  Thus,  the  expression  for 


50 


the  relative  error  in  becomes 

,  r  i^T  (^f (21) 


where  the  uncertainties  in  C,  I,  and  R  are  assumed  to  be 
independent.  For  the  determination,  two  types  of  range 

error  were  assxamed  to  exist;  one  because  of  ship  drift  during 
the  time  data  was  being  collected,  and  the  second  because  of 
surveying  uncertainties.  Therefore,  becomes 


- 

7^ 


(22) 


where  (^)  is  the  relative  error  introduced  by  assxjming  a 
20  meter  ship  drift,  and  the  (^~^qqq)  term  arises  from 
assuming  a  300  meter  error  at  6000  meters  and  a  lesser 
error  as  range  decreases.  Furthermore,  it  is  assumed  that 
=  0.05,  and  a  value  of  5.0  x  10  ^  is  utilized 

c  T 

for  . 

Then,  utilizing  equation  (21),  the  above  assumed  values, 
and  ranges  from  1000  to  6000  meters  in  1000  meter  increments, 
a  plot  of  the  relative  error  versus  range  was  constructed 
(Fig.  14). 

In  addition  to  the  above  investigation,  the  least  squares 
calculations  used  in  the  curve  fitting  process  were  used  to 
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determine  the  variance  in  the  slope  of  the  plotted  data. 
The  Standard  Deviations  of  the  data  for  28  April  and  29 
April  were  found  through  separate  calculations;  the 
variance  of  each  day's  results  were  then  computed. 

This  procedure  indicated  a  relative  error  in yx,  of 

167o. 
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VI.  SUGGESTIONS  FOR  FURTHER  INVESTIGATION 


The  following  paragraphs  are  intended  to  present  several 
ideas  and  modifications  that,  if  implemented,  would  allow  a 
more  accurate  determination  of  the  extinction  coefficient 
under  the  conditions  experienced  during  this  type  of  experi¬ 
ment.  Furthermore,  it  would  be  possible  to  investigate 
the  propagation  of  laser  produced  radiation  at  other  wave¬ 
lengths  . 

A  new  scanning  device,  recently  procured,  should  serve 
to  stabilize  the  "Fan-Scan"  output  under  varying  line  volt¬ 
age  conditions . 

A  second  demodulator  unit,  presently  under  construction, 
will  provide  an  input  to  a  digital  voltmeter,  thus  facili¬ 
tating  a  more  accurate  and  rapid  determination  of  the  laser 
source  intensity.  Additionally,  a  method  of  monitoring 
the  signal  detection  equipment  performance  would  prove 
valuable. 

The  use  of  a  Laser  Rangefinder  should  be  considered. 

This  would  provide  more  accurate  and  rapid  source-to- 
detector  range  determination. 

Ship  to  ship,  or  ship  to  some  other  waterborne  platform, 
operations  should  be  considered.  Since  the  shoreline,  across 
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which  the  propagation  path  passes  in  the  present  experi¬ 
mental  setup,  is  rock  strewn  and  relatively  unprotected, 
waves  breaking  on  these  rocks  produce  a  significant  amount 
of  airborne  water  particles  and  vapor.  This  being  directly 
in  front  of  the  source  can  only  act  to  produce  a  fallacious 
result  in  the  determination  of  extinction.  An  experiment 
conducted  entirely  over  water  would  eliminate  this  problem. 

Further  attempts  should  be  made  to  measure  0.6328  micro¬ 
meter  radiation  extinction  in  fog,  mist,  and  rain.  Also, 
an  effort  should  be  made  to  measure  the  extinction  coeffi¬ 
cient  of  other  wavelength  laser  radiation;  specifically, 

CO2  produced,  10.6  micrometer  radiation. 
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APPENDIX  A 


SPHERICAL  WAVE  SOLUTION  TO  THE  WAVE  EQUATION 

The  following  development  will  show  that  a  laser,  which 
produces  pure  monochromatic  radiation,  can  be  considered  to 
be  a  point  source  bf  spherical  waves  as  long  as  the  region 
of  observation  of  these  waves  is  restricted  to  a  region 
near  the  beam  axis  and  at  a  range  from  the  source  much 
greater  than  the  distance  off  the  axis. 

Consider  now  a  laser  source  in  which  the  intensity 
distribution  at  planes  normal  to  the  propagation  is 
Gaussian.  Starting  with  the  Maxwell  Equations  in  a  homo¬ 
geneous  charge-free  medium,  so  that  V-  E=0  ,  the 
characteristics  of  the  Gaussian  Beam  are  developed  /Ref.  27* 

Then  7  X  ^  ^  ^  J  y  (A-1) 


Where 


T  -  o—  £ , 


Also 


Vxe  =  J// 


(A-2) 


c)€ 


Taking  the  curl  of  (A-2)  and  substituting  (A-1)  for 
YkH  results  in 
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and 


(A-4) 

(A-5) 

where  Y  ^  ^  x  S  =  T  (  V-  ^ )  -  Y  V*  ^ 
and  V  *  ^  -  O  , 


7  <  (Vx  t)  --  ^  i-^  j§-  > 

v"e  =°  . 


Assuming  that  the  field  quantities  vary  as  ^  = 

7^<  S/jz)  -  that  is  pure  monochromatic  radiation 

such  as  from  a  laser  -  the  Wave  Equation  (A-5)  becomes 

Y^£  -hk^Cr)^  =  o  J  (A-6) 

where  ^  6  [  I  -  J 

and  (A-6)  allows  for  the  possible  dependence  of  k  on 
position  r.  At  this  point  the  development  is  limited  to 
the  case  in  which  k^(r)  is  given  by 

k^  (  r  )  =  k^  -  kk2r^  ,  (A-7) 

where  k2  is  some  constant  and  k  =  ,  where  A  is  the 

wavelength  of  a  TEM  (Transverse  Electromagnetic)  plane 
wave  propagating  in  the  medium.  Considering  the  propaga¬ 
tion  of  a  nearly  plane  wave  in  which  the  flow  of  energy  is 
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predominantly  along  the  S-  axis  and  assuming  a  solution  to 

(A-6)  whose  transverse  dependence  is  on  r  =  +  y^  only. 

2 

then  Y  in  (A-6)  can  be  replaced  by 


T"’-.  9;^  +  4—  =  _  +  jL  j' 


taking  E  as 

^  -e 

equation  (A-6)  becomes 


— 


(A-8) 


(A-9) 


-Mtz 


and  operating  on  (A-10)  with  (A-8) 


where 


^'.j4 


,  and  the  assumption  is  made 


that  the  transverse  variation  is  slow  enough  that 


f 

Taking  ^  in  the  form  of 


f  =  [7c^)  r'^J }  ^ 


(A-12) 


substituting  into  (A-11) ,  and  using  (A-8)  gives 


-  kr^q’-2kP'  -  kk- 


r'^-  O. 


(A-13) 


58 


If  (A-13)  is  to  hold  for  all  r,  the  coefficients  of  the 
different  powers  of  r  must  be  equal  to  zero.  Thus 


+  kq''  +  kk2  =  0,  (A-14) 

and  .  (A-15) 

^  k 

Now  according  to  (A-7) ,  if  the  medium  is  homogeneous,  k2 
can  be  set  equal  to  zero  and  (A-14)  becomes 

q2''+  kq'  =0.  (A-16) 


Now  if  the  function  S(z)  is  introduced  by  the  relation 


Q  = 


(A-17) 


It  follows  from  (A-16)  that  S”  =  0  and  s'  =  a  also 
S=  az  +  b  so  that  from  (A-17) 

q(Z)  =  k  .  (A-18) 

Following  Kogelnik,  /Ref.  47,  it  will  be  advantageous  to 
define  a  parameter  q  (z)  where 


St 


Now  (A-18)  may  be  rewritten  as 


(A-19) 


Equation  (A-15)  may  now  be  rewritten  as 


(A-20) 


(A-21) 


59 


and 


>  -^fr)  ) 

where  the  arbitrary  constant  of  integration  is  chosen  as 
zero,  which  amounts  to  fixing  the  phase  of  the  final  solu¬ 
tion  at  zero.  Substituting  (A-20)  and  (A-22)  into  (A-12) 

^  Vj .  (a-23) 

Taking  the  arbitrary  constant  of  integration  to  be  purely 
imaginary  which  leads  to  waves  whose  energy  density  is  con¬ 
fined  to  a  region  near  the  Z  axis,  and  re-expressing  it  as 
a  new  constant  cJo 

f.  --  ,  (A-24) 

where 

--  .  (A-25) 

'v 

Using  (A-24)  and  the  fact  that 

the  first  term  of  (A-23)  becomes 


J'i  ‘ ^ 'fr '  jp~h 


{k-ie) 
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Substituting  (A-24)  into  the  second  term  of  (A-23)  and 
separating  into  real  and  imaginary  parts 


-r 


(A-27) 


Defining  the  following  parameters 


-  (.)  •  «•// .  (it,Tl , 

(A-28) 

Tf  .  ,1,  ^5-77  , 

(A-29) 

(4Aj>)  , 

(A-30) 

Combining  (A-26)  and  (A-27)  into  (A-23) ,  and  utilizing  the 

fact  that  E(;^^^*)  =  €.  results  in 

i)  -  (k€.-4>)  -  r  * 

jr)J 

Consider  now  the  form  of  a  spherical  wave  emitted  by  a 
point  radiator  placed  at  Z  =  0  /Ref.  37.  It  is  given  by 

=<  -^■e.7c/>  (-^k 


«<  :k  -J,k  ^  ) 


(A-32) 


?  9  2 

provided  that  X  +  Y  <■<  Z  ,  the  paraxial  limitation,  is 
observed. 

The  last  form  of  equation  (A-32)  follows  where  Z  is 

equal  to  R,  the  radius  of  curvature  of  the  spherical  wave. 

?  2  2 

Since  r"^  =  X  +  y  ,  equation  (A-32)  becomes 

-  r  3  •  (A-33) 

Now,  comparing  the  form  of  equation  (A-31)  to  that  of 
(A-33) ,  it  can  be  seen  that  (A-31)  is  in  the  form  of  a 
spherical  wave,  provided  the  observer  is  in  the  paraxial 
region  at  a  distance  much  greater  than  the  off-axis 
distance,  and  also  provided  that  ■  >;>  1.  . 
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APPENDIX  B 


PROCEDURE  FOR  DETERMINING  AVERAGE  VOLTAGE,  tT  , 
FROM  LOG  VOLTAGE  DISTRIBUTION 


The  following  development  will  serve  to  illustrate  the 
validity  of  taking  the  peak  value  of  the  nearly  Gaussian 
Log  voltage  distribution,  as  modified  by  the  standard  devia¬ 
tion  of  the  best  fit  Gaussian  curve,  to  be  a  measure  of  the 
average  observed  voltage  . 


Let 


where 


vz 


(B-1) 


channel  number  in  multichannel  analyzer, 
voltage  into  the  log  voltage  circuit. 


and 


/(  =  0^2 '2/, 


(  XT,  and  Xtx,  known  calibration 
voltages) , 


7/'  =  -  '0~i 


Now,  define 


So  that 


•r--  ^/H 


(B-2) 


(B-3) 
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and 


From  (B-3)  and  (B-4) , 


J-X-- 


Xt  -  Vl  e.  . 


Now 


^  pL'ic)cL'?c 


where 


^TL'^)cL'^  -  ^  p{,'K)cL'X 


Thus,  the  average  voltage  can  be  expressed  as 


—  ^  V~ 

^  d-d-  ) 

or  using  (B-7) , 

'U'  -  ^  V-cl'O'^ 

from  (B-6) , 

V--  -ir  f>  (.■>:) 

Substituting  (B-5)  into  (B-10)  for  'IT  results  in 


(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 


Now  assxame 


^  /f  e. 


rr-X.l^ 

^2  0^- 


(B-ll) 


(B-12) 
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and  replacing  p(x)  in  (B-11)  results  in 


ir  -- 


r  -  L-r-  ^ 


(B-13) 


or 


v-  - 


-C±Ml^Af^l2d 

'iTo!^  \  ^  A  y 

(B-14) 


Expanding  the  exponential: 


(B-15) 


Replacing  the  exponent  in  (B-14)  by  (B-15)  results  in 


ji(nc 


(B-16) 


Since  ^  is  independent  of  %  ,  it  can  be  removed 

from  under  the  integral,  and  recognizing  that  h]€. 


is  equal  to  C,  from  (B-7)  and  (B-12) ,  equation  (B-16) 
becomes 


IK  ^ 


-h  ^ 
r  X 


(B-17) 


In  actual  operation,  p(x)  was  the  distribution  collected  by 
the  multichannel  analyzer.  In  calculating  the  best  Gaussian 
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curve  fitted  to  the  data  distribution,  the  computer 
generated  and  subsequently  V”  .  The  validity 

of  this  result  depends  solely  on  how  well  the  computed 
Gaussian  actually  fits  the  experimental  distribution.  This 
was  indicated  graphically  by  the  curve  drawn  onto  the  data 
(see  Figure  6) . 

Reference  8  describes  the  computer  program  in  detail. 
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APPENDIX  C 


DATA 


SPANAGEL  HALL  BASEMENT 

2  February  1976 

RANGE  (meters) 

IndR^) 

31 

6.89 

62 

7.15 

100 

7.15 

122 

7.31 

150 

7.20 

200 

7.15 

250 

6.95 

275 

6.80 

NAVY  ANNEX 

18  March  1976 

91 

7.23 

91 

7.26 

181 

7.01 

181 

6.83 

181 

6.99 

335 

7.04 

335 

7.06 
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NAVY  ANNEX 

19-20  April  1976 

RANGE 

(meters) 

In(IR^) 

91 

8.19 

91 

8.22 

91 

8.25 

91 

8.15 

206 

8.14 

206 

8.30 

336 

7.99 

336 

8.34 

336 

8.14 

PT.  PINOS 

to  ACANIA 

30  March  1976 

RANGE 

(Meters) 

IndR^) 

1604 

5.42 

2384 

5.69 

3323 

5.19 

3915 

4.74 
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PT.  PINOS  to  ACANIA 

28  April  1976 

RANGE  (meters) 

In(IR^) 

1279 

7.52 

2195 

7.53 

3288 

7.42 

4420 

6.60 

4535 

6.65 

5332 

5.75 

5425 

5.69 

6748 

5.81 

PT.  PINOS  to  ACANIA 

29  April  1976 

RANGE  (meters) 

IndR^) 

839 

7.82 

2120 

7.76 

3251 

7.36 

*4156 

7.35 

5300 

6.12 

6283 

5.89 

■*The  following  is  an  example  of  the  calculations  and 
corrections  applied  to  the  value  of  It  ,  obtained  from 
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the  pulse  height  distribution,  to  arrive  at  the  value  of 
In  (IR^)  tabulated  in  this  Appendix.  The  particular  data 
point  used  in  this  example  was  chosen,  since  all  correc- 

O 

tions  possible  were  applied  in  the  determination  of  In  (IR"^) 


In 


J  2  ?  _  ,  f Average  voltage  from  pulse  height  distribution 

IIR  J  -  In  ^  Pre-amnlifier  Gain  Setting 


31  mv 

Monitor 

Reading 


Pre-amplifier  Gain  Setting 

0.04  .9') 

(Aperture  Radius)^  j 


X 


in  ^ 

,  0.538 

31  0.04 

100 

18.5  ^  4  ^ 

In  1 

>  1557.1244 

] 

7.35 

X  (4156) 
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